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Ab initio real-space Hartree-Fock and correlated approach to optical dielectric 

constants of crystalline insulators 

Priya Sony and Alok Shukla 
Physics Department, Indian Institute of Technology, Powai, Mumbai 400076, INDIA 

In this paper we present an approach aimed at calculating the optical dielectric constant (eoo) 
of crystalline insulators both at the Hartree-Fock, and correlated levels. Our scheme employs a 
real-space methodology, employing Wannier functions as the basic building blocks. The scheme has 
been applied to compute eca for LiF, LiCl, MgO, and Li20. In all the cases, results of correlated 
calculations, based upon systematic many-body corrections beyond Hartree-Fock, exhibit excellent 
agreement with the experiments. To the best of our knowledge, the calculations presented here are 
the first of their kind for bulk solids. 

PACS numbers: 77.22.-d, 71.10.-w, 71.15.-m 
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INTRODUCTION 



It is of considerable interest to be able to compute 
various properties of solids using an ab initio method- 
ology based upon their many-particle wave functionsi. 
Such an approach allows one to initiate the calculations 
at the Hartree-Fock level, and then systematically im- 
prove them by including the electron-correlation effects 
to the desired level of sophisticationi. In our recent 
work, we have developed such a methodology for per- 
forming electronic structure calculations on crystalline 
insulatorsSi^ii^. The approach employs a real-space 
philosophy with Wannier functions (WFs) used as the 
single-particle basia^^^ii. Electron-correlation effects are 
subsequently included by employing the "incremental 
scheme" which is nothing but a Bethe-Goldstone-like 
many-body hierarchy, consisting of virtual excitations 
into localized unoccupied orbitals^^. Early applications 
of this approach consisted of calculations of ground state 
geometries, cohesive energies, and elastic properties of 
crystalline insulators at the Hartree-Fock level^iii, as 
well as at the many-body level^. Recently, we extended 
the approach to perform both Hartree-Fock and many- 
body calculations of the Born effective charges of insula- 
tors using a finite-field approach^. In the present work 
we further extend our approach to compute a very im- 
portant quantity associated with dielectrics, namely the 
optical dielectric constant eoo- Optical dielectric con- 
stant is a quantity to which electron correlation effects 
can make significant contributions. The fact that it is 
difficult to compute this quantity accurately is obvious 
because it is only in the second-order perturbation ex- 
pansion of the total energy that the corrections due to 
the polarizability of the system make their first appear- 
ance. In this paper we demonstrate that our ab initio 
Wannier-function-based methodology can also be used to 
perform accurate calculations of Eqc , by reporting results 
on bulk LiF, LiCl, MgO, and Li20. The theoretical re- 
sults in all the cases are in very good agreement with the 
experiments. Besides exploring the infiuence of electron- 
correlation effects on eoo, we also examine the infiuence 
of the choice of basis functions on this quantity. We con- 



clude that the best computation of electron-correlation 
effects is achieved only if the basis set is of good quality. 
Additionally, wherever possible, we also compare our re- 
sults on eoo to those obtained using the density-functional 
theory (DFT) based approaches. 

Remainder of this paper is organized as follows. In 
section^lwe briefiy discuss the theoretical aspects of our 
approach. Next, in section ITTTI we present and discuss the 
results of our ab initio calculations performed on various 
systems. Finally, in section IIVI we present our conclu- 
sions. 



II. THEORY 

Starting point of our approach is the expression of the 
electronic contribution to the polarization per unit cell 
(dipole moment per unit volume, P^'^-') defined as the 
expectation valueifi 



P^^^=^(*0|Re|*0), 



(1) 



where fl is the volume of the unit cell, qe is the elec- 
tronic charge, N (— > cx)), represents the total number 
of unit cells in the crystal, Rg — J2k=i ^^ ^^ ^^^ many- 
particle position operator for the N^ electrons of the crys- 
tal, and l^o) represents the correlated ground-state wave- 
function of the infinite solid. From the polarization, the 
dielectric susceptibility tensor (per unit cell) of the solid 
iXij) can be computed as 
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where, E is the external electric field, and i and j refer 
to the Cartesian components of the vector quantities in- 
volved. For the cubic crystals considered here, only the 
diagonal components of Xij are nonzero, which we de- 
note as X; henceforth. Finally, the dielectric constant of 
the crystal can be computed from the well know expres- 
sion e = 1 + 47rx. By computing the polarization of the 



crystal for atomic positions corresponding to its ground 
state (i.e. no relaxation in the presence of applied electric 
field) , the dielectric constant obtained will be the optical 
dielectric constant Coo of the crystal. In the present work, 
we avoid the tedious route of calculating first the corre- 
lated many-particle wave function of the system (|*o)), 
followed by the expectation value of the dipole operator 
(cf. Eq. (QJ) and its subsequent derivative. Instead, 
we use the generalized Hellman-Feynman theorem, and 
the finite-field approach of computing dipole expectation 
valuesii. Accordingly, we employ the modified Hamilto- 
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H'(E) - Ho-QeE 



fc=l 



(3) 



where Hg is the usual Born-Oppenheimer Hamiltonian 
for the solid, and E is a user specified external electric 
fieldi^. This leads to the modified Hartree-Fock equa- 
tions in the Wannier representation 

{T + U + 2J-K + XPS- qeE-r)\a) = e„|a), (4) 

where T represents the kinetic-energy operator, U is the 
electron-nucleus potential energy, J is the Coulomb term 
(or the direct term) of the electron-electron repulsion, K 
is the corresponding exchange term, Ps is the projection 
operator constructed from the WFs in the short-range 
environment £ of the reference unit cell, A is a large shift 
parameter, \a) is one of the WFs in the reference unit 
cell, and ea is the corresponding eigenvalue. The projec- 
tion operator Ps, coupled with A, plays the role of the 
localizing potential, leading to self-consistent calculation 
of WFs localized in the reference unit celliSii The nota- 
tions used here are consistent with those in our original 
work aimed at direct determination of WFs of crystalline 
insulators^"-, and are explained there in detail. The nov- 
elty of Eq. Q , as compared to the original Hartree-Fock 
equationsSii, is the presence of the term corresponding to 
the electric field. Thus the WFs \a) above, are obtained 
in the presence of an external electric field, and there- 
fore include the effects of dielectric polarization as also 
enunciated by Nunes and Vanderbilt for tight-binding 
Hamiltonianl^. One has to be careful in choosing the 
values of the electric field E, because very large values 
can destroy the lower-bounds of the spectrum leading to 
catastrophic resultaiS. In the present work the maximum 
magnitude |E| = 0.01 a.u. was used, and it did not cause 
any numerical problems. 

Once the Hartree-Fock equations for the system — with 
and without the external electric field E — have been 
solved, the correlation effects are included by adopting 
the approach described in our earlier papers i^i^i^ For both 
types of many-body calculations (i.e. with and without 
the electric field), the correlation energy per unit cell is 
computed as per "incremental expansion'—. 
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E ^'^■yfe 
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where e^, Ae^, Ae^fe, . . . etc. are respectively the one-, 
two- and three-body, . . . correlation increments obtained 
by considering simultaneous virtual excitations from one, 
two, or three occupied WFs, and i, j, k, . . . label the 
WFs involved?. The total correlated energy /cell Eceii can 
be simply obtained by adding Ecorr to the Hartree-Fock 
energy Ehf of the system, i.e., Eceiii^) = Ehf(Ei) + 
Ecorr (El)- In the previous equation we have included the 
external electric field E expHcitly in the parentheses to 
highlight the fact that all the energies involved in the ex- 
pression depend upon it. As in our earlier works, the cor- 
relation approach used for computing various increments 
was the full configuration-interaction (FCI) method, and 
the incremental expansion was restricted to the two-body 
terma^iM. For all the materials investigated, there are 
four valence WFs (ns and three np functions, n being 
the principal quantum number of the valence band) per 
unit cell. Therefore, the number of one- and two-body 
increments for these materials is quite large. However, by 
using various point-group symmetries, we significantly re- 
duced the total number of increments actually computed. 
Finally, the susceptibility/cell of the system concerned is 
obtained using the relation 
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cell 
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The second derivative of the total energy/cell, with re- 
spect to the external electric field needed in the equa- 
tion above is computed numerically, using the well-known 
central-difference formula. 



III. CALCULATIONS AND RESULTS 

Present calculations were performed using lobe-type 
Gaussian basis functions employed in our earlier works 
as well2i2iiiMi&. For MgO, we used the basis set used 
by McCarthy and Harrison^i. For Li20, LiCl, and MgO 
we went beyond the sp-type basis sets, and augmented 
them with a single d-type function centered on the an- 
ion, with the exponents 0.8 for O in MgO, 0.65 for O 
in Li20, and 0.9 for CI in LiCl. The short range en- 
vironment region £ (cf. Eq.l^l) was taken to include 
up to third-nearest neighbors of the reference unit cell. 
Therefore, the unit cell WFs for all the materials stud- 
ied were described using basis functions centered in the 
cells as far as the third-nearest neighbors of the reference 
cell2*^. For all the compound we used the experimental 
values of the lattice constants which were of 4.19 A for 
MgO, 4.573 A for LiaO, 3.99 A for LiF, and 5.07 A 
for LiCl. For all the systems, the fee geometry was as- 
sumed. In MgO, LiF, and LiCl, the anion and the cation 
locations in the primitive cell were taken to be (0, 0, 0) 
and (a/2, 0, 0), respectively, while in Li20, the anion was 
located at (0,0,0), and cations were located at positions 



(±a/4, ±a/4, ±a/4), where a is the lattice constant. Dur- 
ing the correlated calculations, for LiF and Li20, the Is 
WFs of both the cations and the anions were treated as 
frozen core. In MgO, Is WFs of both the atoms, while 
2s, and 2p WFs located on the Mg were held frozen. In 
case of LiCl, Is WF of Li, and Is, 2s, and 2p WFs of 
CI were held frozen, while 3s, and 3p chlorine WFs were 
correlated. Thus, during the correlated calculations, the 
virtual excitations were considered only from the valence 
WFs of the respective compounds. For numerical cal- 
culations of the susceptibility (see Eq. Q), we used the 
central difference formula for the second derivative, along 
with the values of the external electric field E ranging 
from ±0.001 a.u. to ±0.01, in the x direction. 

Results of our calculations are presented in tables ID 
and|n| Finally, in table lllll our best results are compared 
with the experimental values, and also with those of other 
authors. First, we examine Table ID which presents the 
calculated values of Coo for LiF, Li20, LiCl, and MgO 
obtained at the HF level, and at various levels of correla- 
tion treatment, employing different basis sets. As men- 
tioned earher, we have included correlation corrections 
using a real-space Wannier-function-based approach de- 
veloped earlier by us^i^. In the present case, the in- 
cluded correlation contributions correspond to one- and 
two-body increments confined to the WFs in the refer- 
ence cell. Thus two-body contributions corresponding to 
WFs in the nearest-neighboring cells (and beyond) have 
not been included. Elsewhere, we will report that the 
correlation contributions to Coo arising from two-body in- 
crements decrease rapidly with the increasing distances 
between the WFs involved. Therefore, the inclusion of 
correlation increments reported here will suffice for most 
cases. When we compare the results of our best calcu- 
lations to the experimental ones (cf. table Hill , we find 
that generally the agreement between the two is excellent. 
The maximum disagreement is for MgO, with a relative 
error of 6.8%. For other three compounds the relative 
error is well below 5%. Next we analyze the results of 
these materials one-by-one. Upon comparing the results 
of our Hartree-Fock values of eoo for LiF with the experi- 
mental one, we obtain very good agreement between the 
two with the theoretical value being slightly smaller than 
the experimental one. Next, when we include the corre- 
lation effects, the theoretical value increases a bit, but 
still stays in a excellent agreement with the experimental 
value. For the case of MgO, we note that the results ob- 
tained with the sp basis set are in poor agreement with 
the experiments, the disagreement being «27%. Upon 
including the single d-type exponent, the improvement 
in the results begins at the HF level itself, with the cor- 
relation contributions eventually bringing the results in 
close agreement with the experiments. In tabled similar 
trends are seen for the remaining compounds Li2 and 
LiCl in that the results are generally in poor agreement 
with the experiments if only the sp basis sets are used, 
however, inclusion of d exponent brings them in excel- 
lent agreement with the experiments. Noteworthy point 



is that it is only with the inclusion of d exponents that 
correlation contributions are significant, a point which is 
well-known in the quantum-chemistry community in the 
context of polarizability calculations^. Now the ques- 
tion arises as to why excellent results were obtained for 
LiF just with the sp basis set, while for MgO, Li20, and 
LiCl inclusion of the d exponents was essential. The fact 
that for LiF, excellent results are obtained even without 
the use of d-type functions, is a clear indication of the 
fact that in F~ ion valence electrons are extremely local- 
ized owing to the large electron affinity of fiuorine, and, 
therefore, can be described well with an sp basis set. For 
Li2 we attribute the large infiuence of d-type functions 
to the extremely diffuse nature of the valence electrons 
of O ion. As a matter of fact, 0~~ ion does not exist 
in free form, and is stabilized only because of the crystal 
field. Similarly, Cl~ ion is comparatively more diffuse as 
compared to the F~ ion, thereby requiring the use of the 
d-type functions. The other trend which is clearly visible 
from tabled is that when the calculations are performed 
with a proper basis set, inclusion of the electron correla- 
tion effects brings the results in much closer agreement 
with the experimental values, as compared to the HF 
results. 

Next in table lUl we examine the relative contributions 
different correlation increments to the value of eoo of 
these materials. Inspection of the table reveals that some 
very clear trends are visible regarding various contribu- 
tions. In all the cases it is seen that the inclusion of 
the one-body increments increases the value of the eoo 
with respect to the HF value. As far as the contribution 
of the two-body increments is concerned, tahlelTTI reveals 
that for all the materials except for LiF, this contribution 
is negative when the sp basis set is employed. This, in 
our opinion, is a consequence of an inadequate basis set. 
Once the spd basis sets are employed for these materials, 
the contribution of this increment becomes significantly 
larger in magnitude, and positive in sign, to bring the 
final values of eoo in excellent agreement with the exper- 
iments. 



Finally, in table lllll we compare the results of our cal- 
culations with those of other authors. As far as the DFT- 
based approaches are concerned, one very powerful ap- 
proach aimed at computing eoo is the density-functional 
perturbation theory (DFPT) methodic. But, for the ma- 
terials considered here, we were unable to locate any cal- 
culations of eoo based upon DFPT. However, in the lit- 
erature we found several results on eoo based on other 
DFT-based approaches which we discuss next. Umari 
and PasquarelloiS, presented ah initio molecular dynam- 
ics based calculations of eoo of MgO, performed in the 
presence of a finite electric field. Mei et alM performed 
DFT-based calculations of eoo of LiF and LiCl, using lo- 
calized charge densities. For LiF and several other ma- 
terials, Bernardini and Fiorentiniii reported calculations 
of eoo using a "polarization approach" in which they per- 
formed DFT-based calculations on slabs of the material 
to calculate macroscopic polarization. Then they ob- 



tained the values of Coo by using a relation between the 
interface charge, polarization difference and eoo- To the 
best of our knowledge no prior calculations of eoo of Li2 
exist. When we compare our results on MgO with those 
of Umari and Pasquarelloi^, it is obvious that the two 
results are in excellent agreement with each other. For 
LiF, our results and those of Mei et alM^ agree almost 
perfectly, while the ones of Bernardini and FiorentiniAi 
overestimate the experimental value. For LiCl results of 
Mei et al. overestimate Eoo somewhat as compared to 
the experimental value, while our results are almost in 
perfect agreement with the experiments. 



IV. CONCLUSIONS 



experimental results was obtained. The attractive 
aspect of the approach is that the correlation effects 
are organized as per a Bethe-Goldstone-like many-body 
hierarchy, which makes itself amenable to systematic 
expansion or truncation. The influences of various terms 
in this correlation hierarchy on the value of e^o was 
examined for all the systems, which, in our opinion, 
leads to a deeper understanding of the many-body 
effects. Additionally, the influence of the choice of basis 
set on the computed values of Coo was also explored. 
The results presented in this work suggest that this 
methodology holds promise, and its applications to more 
complex systems such as perovskites will be pursued in 
future. 



In conclusion, we have presented an ah initio wave- 
function-based approach aimed at computing the optical 
dielectric constant of insulators both at the Hartree-Fock 
level, and the correlated level, employing a real-space 
Wannier-function-based formalism. The approach was 
applied to compute the eoo of crystalline LiF, Li20, 
LiCl, and MgO, and very good agreement with the 
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Table II: Contribution of various types of electron-correlation 
effects to too- Both the one-body and the two-body contri- 
butions arise from virtual excitations from the WFs located 
within the reference cell. Heading two-body (O) implies the 
corrections due to simultaneous excitations from two distinct 
Wannier functions located on the anion in the reference unit 
cell. 





(■ao 


System 


one-body 


two-Body(O) 


LiF(sp basis set) 


0.105 


0.015 


MgO (sp basis set) 


0.039 


-0.008 


MgO (spd basis set) 


0.087 


0.338 


Li2 0(sp basis set) 


0.034 


-0.016 


Li20(spd basis set) 


0.239 


0.295 


LiCl (sp basis set) 


0.042 


-0.018 


LiCl (spd basis set) 


0.043 


0.694 



Table III: Comparison of our best values of too of various 
materials with those of other authors, and experiments. 



System 


This work 


Other works 


Experiment 


MgO 


2.76 


2.78" 


2.96"^ 


LiaO 


2.62 


— 


2.68" 


LiF 


2.04 


2.03*, 2.19^= 


1.96^ 


LiCl 


2.83 


2.95'' 


2.79^ 
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